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We discuss properties of the electroweak vacuum as a function of an external magnetic field. The 
interest in these properties arises due to possible existence of the electromagnetically superconduct- 
ing phase of QCD in the background of a sufficiently strong magnetic field. In the electroweak 
theory, a superconducting vacuum exists for a well defined range of magnetic background fields, 
and the interest of the current paper is the behaviour near the upper critical magnetic field. We 
determine this critical field and calculate the vacuum condensates to leading order. 
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I. INTRODUCTION 

In the past few years, interest has been renewed into 
the effect of strong magnetic fields on the vacuum. It 
was pointed out [T] that in the case of strong interac- 
tions, the lightest vector mesons undergo a phase tran- 
sition and condense into an Abrikosov lattice [2j. The 
vector mesons form a new condensed state of vector 
mesons which is similar to the magnetised state of the 
electroweak vacuum. The structure of the condensates 
in the electroweak vacuum was first described by Amb- 
jorn and Olesen [5] more that 20 years ago. Recently it 
was found that the main feature of both - electroweak 
and rho meson - condensed structures is their unusual - 
magnetic field induced - superconductivity 0] [5] . This 
particular kind of superconductivity is instated above a 
critical value of the applied magnetic field and the con- 
densates are reinforced with growing magnetic fields. By 
comparison, in ordinary superconductivity, mediated by 
Cooper pairs, the superconductivity breaks down above 
the critical magnetic field. Both models, for electroweak 
vector bosons and rho mesons, contain a SU(2) multi- 
plet of vector fields, and there's little wonder that the 
results are found to be so similar. The main distinction 
between QCD and the electroweak theory is the scale dif- 
ference: Because the masses of the rho-mesons are about 
100 times smaller than the W-mass, the critical magnetic 
field is 10 4 times weaker. 

A first look at the superconducting regime - involving 
mostly classical calculations - has given some interesting 
results above but very close to the critical magnetic field 
B c \. Working close to the critical field allows to work 
in the limit where the newly created vector condensates 
are very small. The thus found condensates are similar 
for both the superconducting QCD and electroweak vac- 
uum. The main conclusion is that a hexagonal lattice 
of Abrikosov vortices appears as lowest energy state for 
the charged vector condensate. The neutral condensates 
are essentially derived form this hexagonal structure and 
form a more involved pattern. 

For the electroweak vacuum it is known that at very 
high magnetic fields, way above the one for which vector 



boson condensation and superconductivity sets in, an- 
other phase transition occurs [5J. The Higgs condensate 
will vanish at this critical magnetic field, B C 2, similar to 
what one observes at the electroweak phase transition at 
some critical temperature. Indeed, at higher magnetic 
fields, the full SU(2) group is a symmetry of the classical 
vacuum. As we will show, none of the SU(2) gauge fields 
will be present in the vacuum as a condensate. The struc- 
ture of condensates at such enormous magnetic fields are 
relevant for example in baryogenesis scenario's [7J. 

It is the goal of the present paper to give a description 
of the vector condensates in magnetic fields near, this 
second critical magnetic field (region II in Figjl). A de- 
scription was already given in the case of the non-physical 
Bogomolnyi limit [8] , an analysis that can be extended to 
the general case, as we will show. The Higgs condensate 
will be approximately vanishing while the SU(2) group is 
still (softly) broken. This allows for condensates of chro- 
momagnetic vortices while they are absent in the pure 
Yang-Mills case. We give an expression for these conden- 
sates similar (and in some cases equivalent) to the one 
found by Olesen by solving a Liouville equation in 
terms of elliptic Weierstrass functions. 

We start by giving a description of the electroweak the- 
ory in function of the relevant degrees of freedom, those 
transversal to the applied magnetic field. From this we 
find non-trivial solutions to the pure Yang-Mills equa- 
tions, describing approximately the vector condensates 
near B c2 , the region II of Fig[T] Next we explain the 
phase transition in detail and argue that, classically, all 
condensates should become trivial in the symmetry re- 
stored phase. 



II. THE ELECTROWEAK MODEL IN A 
MAGNETIC FIELD 

The bosonic electroweak Lagrangian in its simplest 
form separates SU(2) gauge fields, the Higgs field and 
the U(l) hypercharge fields. 
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FIG. 1. The electroweak vacuum at zero temperature as 
a function of a constant magnetic background. While the 
condensates in the superconducting vacuum obey non-linear 
equations of motion, accurate approximations can be found 
for small vector condensates (region I) and small Higgs con- 
densates (region II). 



While the latter are conventionally denoted by B^, we 
have adopted the notation instead to avoid confusion 
with the magnetic field B. 

A more suitable form is found when defining the elec- 
trically charged fields W£, for which we have to fix 
the direction of the Higgs vacuum expectation value: 
$ = (0, <fi). One finds then: 
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The magnetic background field will be assumed to be 
directed in the third spatial direction and we will refer 
to the (l,2)-plane as the transversal plane. It is clear 
that only polarisations in this plane couple to the back- 
ground, and may form magnetically induced condensates. 
Restricting ourselves to transversal effects, it is more con- 
venient to work with the energy density instead of the La- 
grangian formalism and to use the complex coordinates 
z = x\ + iX2 and z = x± — ix<i and likewise for the fields: 

W = W±+ 1W2, Z = Z X - iZ 2 , S = 2>i - iD 2 . . • (9) 

The vacuum equations of motion are then fully deter- 
mined by the energy density as a function of transversal 



degrees of freedom. 
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An important observation in the above potential ( 10 ) 
is that only the second term depends on the phases of 
W and W, while the rest of the potential reaches its 
minimum at 



W = 0, 



(11) 



for non a non-vanishing Higgs condensate. While it is not 
certain that one can minimise the full potential with (111, 



this seems a reasonable approximation. It is known that 
for small W condensates, B > B c i, solutions assuming 
(11) exist and we will show that vacuum solutions for 
small Higgs condensates, B < B C 2, exist with W = 
and with the phase of W such that 

SW = 0. (12) 



Therefore we can eliminate W from ( 10 1 entirely, and 



with some manipulations one can rewrite it as follows: 



f = H (^2-^-^M|-sin^|W| 2 
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where the left out boundary terms are: 
_ dd\c/)\ 2 - d(<pD<t>) ~ d{4>D$) 



2 cos 9Fi2 — sin 9Z\ 2 1 cos 4 9 . 
+ g-^0 2^9^ (U) 



The special case for which the coefficient of the last term 
in ( 13 1 vanishes, gives an easy solution for all values of 



the magnetic background field, setting 
f i B = M H = M z . 



(15) 



This Bogomolnyi limit was treated in detail in [8]. 

It is clear that for small \<f)\, which case is realised in 
the vacuum for a magnetic background field B < B c2 , 
the fourth order term is irrelevant and the conclusions 
about the phase transition at B c2 are independent of the 
Bogomolnyi limit and valid for physical values of the pa- 
rameters. 
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III. THE VACUUM WITH AN INFINITESIMAL 
HIGGS CONDENSATE 

Keeping Fig[l] in mind, we will be considering the re- 
gion II close to B C 2 on the magnetic axis, where we can 
define a small order parameter: 

c 2 = (H 2 )- (16) 

We will consider only the condensates to leading order 
in e. It is clear that this does not imply the fourth order 
terms in the Higgs field, such that we are left with a sum 
of squares in (13). The first two terms can be minimised 



independently, giving linear equations of motion in the 
vacuum: 



sin0F 12 + cos6»Z 12 = ^\W\ 2 + ^|0| 2 , 



(17) 



cos 6F 12 -sm6Z 12 = C °^| + cot 0~ 10| 2 . (18) 



What is left is a Ginzburg-Landau like model depends on 
the Z-flux and Higgs field only: 

I, |2 , 2 cos6>F 12 -sm9Z 12 
101 "^W (19) 



\m 2 +\(A-M 2 H 



where one of the boundary terms is now explicitly writ- 
ten. The decoupling of the equations of motion into lin- 
ear ones by introducing [is only gives a nonzero leading 
contribution to the Higgs condensate if we take 



Ml = M%. 



We are left then with the equation 

.9 



(20) 



Dcj) 



d_ 

dz 



(-2iZ) = 0. (21) 
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As Z12 is a topological quantity in the presence of a Higgs 
condensate, boundary conditions are independent of lo- 
cal fluctuations. If we solve the above equation for a 
homogeneous Z12 = (Z12), for which we have locally 
—2iZ = z(Z\2), we see that physical, bounded solutions 
are possible only for 



(Z 12 ) < 0. 



(22) 



Now we see that minimising the energy by means of the 
Z-flux is entirely determined by the boundary term in 
( 14 ) and we find 



(Z12) = 0. 



From ( 18 1 , ( 20 1 and ( 23 1 we can now identify the exact 



value of the second critical magnetic field by taking \(f>\ 
0: 



eB c2 = Ms = M H . 



(24) 



This coincides with the result of [5] in the Bogomolnyi 
limit. 



Moreover, one can rewrite Eq.(17) at B C 2 
W? 2 = w\ 2 - 9 -\W\ 2 = 



(25) 



because 4> vanishes there. This gives us a clear picture 
of what happens in a compact region of space where the 
magnetic field exceeds B C 2, Fig(2] At the boundary of 





FIG. 2. Squeezed magnetic flux lines, exceeding the critical 
value B C 2 in a compact region. In the region II, only a very 
small Higgs condensate is present, vanishing at the critical 
magnetic field. In between the full lines, the flux propagates 
as hyperfiux only, and no condensates are present. 



such a region the pure Yang-Mills vacuum equations are 
exactly satisfied according to ( 25 ) , and inside -in absence 
of a Higgs condensate- the equations of motion are SU(2) 
symmetric. Therefore one finds no more vector conden- 



sates in this region at the classical levelQ The magnetic 
flux does therefore propagate as (conserved) U(l) hyper- 
charge flux: 



g'X 



12 



-X 12 = Mfj. 

cos 6 



(26) 



IV. PURE SU(2) VACUUM EQUATIONS 

We will devote the remainder of the paper to deter- 
mining the leading order contribution to the vector con- 
densates present for magnetic fields just below B c2 . As 
we have shown, ( 12 ) and (25 ), this amounts to solving the 



(23) classical vacuum equations for a pure SU(2) theory, 



The presence of a nonzero net flux on the contrary 
would have corresponded to Z-strings [JO] with zeroes 
in the Higgs condensate, organised in an Abrikosov lat- 
tice structure. The absence of Z-strings near B c \ was 
already investigated in [IT] . Now it becomes clear that 
also higher magnetic fields do not spoil this property. 



£)W = 0, gwl 2 = y -\W\ 2 , 



(27) 



1 Quantum effects invalidate this statement. |12| 
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but with SU(2) broken by the small Higgs condensate. 
We can emulate this explicit symmetry breaking by in- 
troducing a boundary condition 



12, 



= eB ext ^ 0. 



(28) 



This problem was first treated in |15) using an Abrikosov- 
like approach, which should be valid for small condensate 
and linearised equations. A more accurate approach was 
given in [9], where solutions for rectangular lattices in 
terms of Weierstrass elliptic functions are obtained. Our 
approach is similar to this, but does cover the case of a 
hexagonal vortex lattice which is expected to be the true 
ground state. 



General solution 



We rewrite (27) as 



BliiW = -igW 3 , 
dlnW* = igW 3 , 
ig(BW 3 - BW 3 ) = -g 2 \W\ 2 , 



from which we derive: 



Aln|W/| 2 + g 2 \W\ 2 = 0, 
A axg(W) + | (BW 3 + BW 3 ) = 0. 



(29) 
(30) 
(31) 



(32) 
(33) 



Choosing the gauge 8 ■ W 3 = BW 3 + BW 3 = makes 
the second term of the second equation disappear and 
aig{W) should be a harmonic function. Naively, one 
would say that the residual gauge symmetry would al- 
low for any such function, but one must keep in mind 
that the non trivial topology of the vortices gives further 
restrictions. Taking that into account will allow us to 
determine the W condensate up to "small" gauge trans- 
formations only. 



Now we turn to ( 32 ) , which is called a Liouville equa- 



tion and has solutions of the form 113 



g\w\ 



2V21 



j/(z)t ■ M ■ y(z) ' 



(34) 



while ( 29 ) can be written as 



dz 



\n\W\ = 



-i(- wg (W) + y -W 3 ). 



(37) 



Let us discuss both terms on the right hand side in the 
gauge B-W 3 = dW 3 + 8W 3 = or, because of ((33}, with 



arg(W^) a sum of a holomorphic and anti-holomorphic 
function. 

The singular part of the argument of W with anti- 
vorticity n will in the vicinity of a single vortex be 



arg(VF) « — hi-, 



(38) 



in a gauge where W 3 is non-singular. We choose n in this 
way because we know from previous work [5] that the 
negatively charged W condensates have negative wind- 
ing, while the positively charged W' have positive vor- 
ticity. It will become clear that also in the current ap- 
proach only for this ansatz solutions can be found. We 
thus find that the singular part of the right hand side of 
(371 is given by: 



d ft 

— In \ W\ — - — h finite part, 
dz 2z 



(39) 



such that we find from (36) for the meromorphic E(z): 
In + n 2 residu 



1 



E{z) 



4z 2 



h{z) 



(40) 



with h{z) some undetermined holomorphic function. 
Before proceeding to the solution of (32), we can 



present a slightly different approach to obtaining the 
form (34). As it is a solution to pure gauge equations 



of motion, the solution should be locally pure gauge and 
therefore have the form 



g W *T a = 




= iU- 



(41) 



for some U £ SU{2). Like before we make the restriction 
of this problem to the transversal plane with W = 



0, and use the conditions ( 29p0 ) 



Here M is a hermitian, positive definite 2-by-2 matrix 
with det(M) = 1 and y = (2/1,2/2) are solutions of the 
equation 



y" + l -E{z)y = 0, 



(35) 



with constant Wronskian w and some function E(z), in 
which the boundary conditions are encoded. 
One can prove 



l -E{z) = ^- 2 \n\W\-{j- z \n\W\) 2 , 



(36) 



U~ l dU 



U- 1 BU = 



■^W BlnVW* 



-ainv^ 



din VW 




(42) 
(43) 



The solution of the above system of equations for U 
gives four independent integration 'constants' yi{z) 1 i — 
L...4: 



U 



Wyl(z) 
Wy* 2 (z) 




(44) 
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for which one essentially uses the triangular character of 
the coefficient matrix and the diagonal elements. Unitar- 
ity gives us then 



\yi\ 2 



\m\ 2 



1 

WY 



(45) 



With the relative signs related to the freedom in the 
choice of the square root branch. Now we can calculate 
with w(z) = yi(z)y' 2 (z) - y' 1 (z)y 2 (z): 



U^dU = 



-dhiVW* 







±2Ww*{z) dlnVW"' 




n- 1 an = ( 9hlVw ±2W * W ( Z ) 
-dlnVW 



(46) 
(47) 



and find the constant Wronskian condition 



tirely consistent with the solution (34 1. The yi can be 



found from the yi by diagonalising the matrix M and 
properly normalising and choosing the phase factor to 
match Eq. fl42|43| ). We then have a 577(2) gauge trans- 
formation that locally gauges away the Liouville solution 
Pi. 



B. Lattice solutions 

Now we will try to put the solutions of the Liouville 
equation (34) on a 2 dimensional lattice, by demanding 



double periodicity in the transversal plane such that each 
lattice cell contains a single vortex. From ( 36 1 we find 



that then the E(z) are elliptic functions, with the inter- 
esting consequence that the 'residu' in ( 40 1 must vanish 



for such functions and the holomorphic part must be a 
constant h. It is not hard to to see from (36(37) that 



such a constant can be gauged away by local U(T) gauge 
transformations. Therefore we can state that a vortex 
around \z\ = in a vortex lattice should correspond to: 



1 2n + r. 

-kjiz) = — 

2 W Az 2 



0. 



(48) 



Now the only elliptic function of the form ( 40 1 for every 
vortex in the lattice, with only one singular point per 
lattice cell is proportional to the Weierstrass ^-function: 

f 



T(z) 



lattice 



(49) 



Also no analytic elliptic functions exist, so Viz) is indeed 
unique for our purpose and we can thus rewrite (|40|) as: 



1 , , 2n + n 2 „ , , 
- 2 E{z) = j—V(z)> 



(50) 



We will construct lattices of vortices of unit winding n = 
1 and consider only that case from now on. 

It is worth noting that the solutions constructed in [5] 
over a square lattice correspond to, in our notations: 



(51) 



and with some algebra one finds for these functions 



y'liz) = 1 3 ( V"(z/2) 
Vi (z) 8\2\V'(z/2) 



V"'(z/2) 
V'{z/2) 



;V(z), 



confirming ( 50 ) with n = 1 



(52) 



We will now try to solve ( 34 ) on a hexagonal lattice 
V n {z) = 4V 3 (z) - \~ 6 



V"{z) = 6V 2 (z), 
assuming the y%(z) functions of V(z): 

Vi(z) - Y l (4V 3 (z)X 6 ). 
This indeed simplifies the equation p5| to 
[X(l - X)Y!' + (c - X{1 + a + b))Y! - abY,] V = 0, 

with 

X(z) = 4A 6 7? 3 (z), 

and with 



a + b = -, ab 
6 



2n + n 2 



144 



2 
3' 



(53) 
(54) 

(55) 
(56) 

(57) 



which is the hypergeometric differential equation. It is 
known that the equation ( 55 1 has 3 singular points X — 



0, 1 and oo, corresponding to the dual lattice points, edge- 
centers and lattice center points, respectively. 
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FIG. 3. The singular points of the yi(z) both on the full 
hexagonal lattice (left) and restricted to just one lattice 
cell (right): The vortex-centers (white circles) correspond to 
X(z) — oo, X(z) = gives dual lattice points (blue triangles) 
and X[z) = 1 mid-points between centers (green squares). 

We will pick two independent solutions such that the 
first is regular (in the sense that it has a well defined 
derivative) at X = and the second at X = 1 H4l: 



Y 1 (X)= 2 F 1 (a,b;-;X), (58) 
Y 2 {X) = X 2 F 1 (a,b;^;l-X), (59) 

with 2F1 the Gaussian hypergeometric function and w 
such that 



\Y 2 \X)Y 1 (X)-Yl(X)Y 2 (X)\ = \w\\X\ 



(X/4)- 2 / 3 

12y/X - 1 



(60) 



G 



which defines the Wronskian for yi and y 2 . To see this, 
note that from (71) and (56 1 it follows that 



dX 

dz 



(z) = 12\- 1 (X(z)/4) 2 / 3 ^/X(z) - 1. (61) 



One finds for the above choice of F's: 
\w\ « 0.25 . 



(62) 



These choices of Yi (X) and thus of yi (z) in p4| give a 
continuous solution over the lattice for any valid choice 
of the matrix M, but extra conditions are necessary to 
get a smooth solution. To cancel the singular derivative 
Y{(1) in the derivative of \W\ we need to ask 



(1, 0) • M ■ Y(l) = 0, 



(63) 



a condition on the first row of M . Similarly the singular 
r 2 '(0) gives 



(0, 1) • M ■ Y(0) = 0, 



(64) 



a condition on the second row. As M is also required to 
be hermitic, to have det(Af) = 1 and be positive definite, 
we have fixed all its degrees of freedom. The result is 
then: 



M = 



r 2 (i)y 2 (o) -Fi(i)r 2 (o) 

-Y 2 (0)yi(l) Yi(0)Yi(l) 

y/WWPWSW) - P)n(i)F 



(65) 



Finally, the scale A is determined by the flux quanti- 
zation condition 



2tt = 



gw\ 2 dx 2 = A(wl 2 ), 



(66) 



with A the surface of one lattice cell given by A = 2^ 
if u>\ is the half-period in the real direction. For the 
Weierstrass functions given by ( 71 ) one can show [H] 



Ar 3 (1/3) = Attuji and so A is given by 



A^r 3 (l/3) 
2vr 



2tt 



9(^12) 




As for the phase, we need some function cr(z) such that 



d 2 in d 2 a*(z) 

— arg(W) = -irTr^ln^- 
oz z 2 az z <j(z) 



-V{z), (68) 



which is by definition the Weierstrass er-function and then 
arg(VF) ee -arg(cr). 

In conclusion, the hexagonal lattice of chromomagnetic 
anti- vortices with unit winding is: 

„ r , , 2721^1 <J*(z) 

gW{z,z)= v ' 1 / > 69 

|yt(z) • M -y{z)\ \a(z)\ 

with yi{z) = y(47 ,3 (z)A 6 ), V and a the Weierstrass T 5 
and a functions, and with the constant M and w given 
by p5| and ((62). 



It is remarkable that the profile found in this way is 
very similar to the one obtained by a more straightfor- 
ward linearised approach with the ansatz 



4 11 /(*), 




(70) 



FIG. 4. The absolute value of the W- condensate as function 
of coordinates on the transversal plane. The dark-to-light fill 
indicates the rising phase, marking the branch points. The 
scale is determined by the applied magnetic field, chosen here 
as g(w\ 2 ) = (5/4) 4 M^. 

Our main goal was to obtain a hexagonal structure, 
but one can prove that for Weierstrass functions obeying 

V' 2 (z) = 4V 3 (z) - \- 4 V(z), V"{z) = 6V 2 (z) - A~ 4 , 

(71) 

and corresponding to a square lattice, the same hyperge- 
ometric equation ( 55 ) as before is found with 



x = Av{zyx t , 



-1/8, 



3/8, c = 3/4, (72) 



instead of (56 1 and ( 57 ) . Repeating the rest of the con- 



struction leads to a square lattice solution identical to 
[5] , but here written in more complicated algebraic form. 



C. Stability 

We have constructed the lattice solutions in the as- 
sumption that SU(2) breaking by a tiny Higgs allowed 
us to fix (w 3 2 ) to a non-zero value, given by the mag- 
netic flux. This situation changes for vanishing Higgs 
condensates, beyond £? c2 . 
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Olesen argued that this lattice, as a solution of the 
pure gauge equation (27), could be stable in a SU(2) in- 



variant theory only if each vortex would have non-trivial 
topology in the SU(2) group. As m(SU(2)) = Z 2 , this 
corresponds to the existence of Alice strings [16 . This 
question can be answered by calculating the sign of the 
Wilson loop over one lattice cell, on^Jcan show that the 
magnetic flux and the condensate contribute compensat- 
ing signs to the Wilson loop and therefore the vortex does 
not carry topological charge. The vortices can therefore 
be gauged away once the symmetry is restored. 

We can recover Olesen's result if we look at the form 
of U in (44). There we actually constructed the gauge 



transformation that makes the condensate and the mag- 
netic flux cancel, which is part of the SU(2) symmetry 
group if no singularities appear. One can easily check 
that, given the expressions for the y or y and for the W 
functions, singularities in (44) could only occur at the 



zeroes of W. However, the most singular of the j/j will 
behave like (W)~? there and thus the singularities can- 
cel against the y/W factor. Therefore, no vortices and no 
chromomagnetic flux appear in the symmetry restored 
phase, and all flux is carried by the hyperflux X\2. 

Also below B C 2 a degeneration of the vacuum lattice 
structure seems to appear in our approximation. The 
fact that multiple lattice structure solve the pure SU (2) 
equations is not surprising. It is known that in the Bo- 
gomolnyi limit of the electroweak model a similar de- 
generacy appears. In that case the fourth order terms 
depending on the Abrikosov parameter [2] for the vortex 
lattice are absent and the energy is a sum of squares. 
The pure SU(2) energy density has the same property. 
Rather than a fixed lattice structure, vortices are free as 
long as they have fixed density given by the flux quanti- 
sation condition) 66 ) . The vacuum is then a vortex liquid. 



While in the presence of a small symmetry breaking 
Higgs condensate the obtained solutions are still good ap- 
proximations, the liquid property will be lost. All previ- 
ous calculations of lattices in the superconducting phase 



between B c \ and B C 2, shown in Fig (T] indicate that then 
the hexagonal structure should be energetically more ad- 
vantageous. Therefore, our new solution, while trivial in 
the symmetry restored phase, is the leading contribution 
in the presence of a small Higgs condensate. 



V. CONCLUSIONS 

We have taken a closer look to the phase transition 
that is reached when increasing the magnetic background 
even further in the superconducting electroweak vacuum, 
as depicted in FigjT] It is already known that the su- 
perconducting phase transition at B — B c \ changes the 
vacuum. This superconducting phase exists in the region 
of magnetic field strength B cl < B < B c2 , and has the 
form of an Abrikosov lattice with a hexagonal structure, 
in which W, Z and Higgs condensates are organised. B C 2 
is reached when the order parameter of the Higgs mech- 
anism, the expectation value of the Higgs field, vanishes. 

We have described the second phase transition at B = 
B C 2 using roughly the same approach as when describ- 
ing the first transition: By considering the condensates 
on one end of the phase transition to be perturbations 
of those on the other end. While a similar approach was 
straightforward when the 'ground state' was the ordinary 
Higgs vacuum, up to now the exact vacuum in the chro- 
momagnetic vortex state was known for a square lattice 
only [5]. 

Furthermore we were able to confirm that qualitatively 
the phase transition proceeds outside the Bogomolnyi 
limit in the same way as was found in [5], and that the 



critical field is given by (24). 



The author thanks Maxim Chernodub, Henri Ver- 
schcldc and Poul Olesen for discussions and suggestions. 
A first version of this paper was submitted unaware of 
[9], many of the points made already appeared there in 
a slightly different form. However, it seems that the ap- 
proaches are sufficiently different to make this work an 
interesting addition. 
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